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Temporal  Attributes  of  the  Ambient  Seismo- 
Acoustic  Environment:  LaJunta,  Colorado 


1,  INTRODUCTION 


In  June  1986,  the  Air  Force  Geophysics  Laboratory  established  a  seml-permanent  selsmo- 
acoustlc  monitoring  array  near  La  Junta,  Colorado,  This  site,  located  under  a  SAC  low-level 
training  route,  was  used  to  Investigate  the  potential  value  of  selsmo-acoustlc  techniques  to 
detect  and  track  low-flying  eilrcraft.  As  a  foundation  for  this  work  It  was  necessary  to 
establish  a  base  measure  of  the  seismic  and  atmospheric  pressure  background  levels  at  the  site 
and  to  define  any  quantitative  relationship  between  the  seismic  and  pressure  noise.  The 
results  provided  In  this  report  are  based  on  the  environment  as  measured  during  the  period  9 
through  29  October  1986. 


2.  RESULTS 


It  Is  the  purpose  of  this  paper  to  establish  the  attributes  of  the  base  seismic  and 
atmospheric  pressure  noise  for  a  surface  array  at  La  Junta.  Colorado,  for  periods  free  of 
conspicuous  events.  Further.  It  was  desired  to  establish  the  level  of  correlation  between 
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surface  pressure  and  seismic  noise  at  this  site.  Due  to  the  nature  of  selsmo-acoustlc  coupling, 
the  main  Interest  is  in  the  relationships  between  vertical  seismic  noise  and  atmospheric 
pressure.  General  statements  of  the  findings  of  this  study  can  be  given  as; 

(1)  The  seismic  and  pressure  noise  levels  observed  at  the  La  Junta  site  in  the  band  1  to  30 
Hz  are  significantly  above  measurement  system  noise; 

(2)  The  seismic  and  pressure  noise  environments  are  neither  Gaussian  nor  stationary  and 
show  distinct  temporal  variations  probably  associated  with  meteorological  conditions: 

(3)  The  maximum  envelope  levels  for  noise  at  this  site,  taken  over  1.01 -sec  windows,  can 
be  adequately  represented  for  the  purposes  of  detection  threshold  estimation  by  Gumbel  Type 
11  extreme-value  distributions,  EXii(u,K).  The  parameters  of  these  distributions  are  [u  =  3.8  x 
10“®,  K  =  2. 56]  and  (u  =  9.1  x  10“^,  K  =  2.29)  for  vertical  and  horizontal  seismlcs.  respectively, 
and  [u  =  0.108,  k  =  2.12]  for  the  atmospheric  pressure  field:  and,  finally, 

(4)  The  seismic  noise  environment  at  the  La  Junta  site  is  only  weakly  correlated  with  the 
total  atmospheric  pressure  signal  in  the  same  band,  and  the  vertical  and  horizontal  seismic 
noise  conditions  are  also  only  weakly  correlated. 


3.  SITE  DESCRIPTION 


3.1  General  Description 

La  Junta.  Colorado  is  located  along  the  Arkansas  River  valley  in  southeastern  Colorado 
and  is  approximately  100  kilometers  east-southeast  of  Pueblo,  Colorado  (Figure  1).  The  site  is 
typical  of  the  high  pleilns  of  eastern  Colorado  with  relatively  shallow  relief  and  semi-arid 
vegetation  consisting  of  sparse  grasses  and  sage  brush.  The  property  surrounding  the  site  is 
used  as  cattle  range  land. 


3.2  Cultural  Setting 

The  towns  of  La  Junta  and  Los  Animas  lie  about  20  km  to  the  west  and  22  km  to  the 
northeast  of  the  site  (Figure  2).  A  multi-lane  highway  trends  east-west  about  13  km  north  of 
the  site  and  runs  along  the  course  of  the  Arkansas  River.  The  highway  is  paralleled  by  an 
active  main  line  of  the  Santa  Fe  Railroad.  A  ramch  service  road  runs  approximately  200  m 
south  of  the  site  and  the  ranch  manager's  home  is  located  about  1.5  km  to  the  southeast. 
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Figure  1 .  Map  Showing  the  General  Location  of  the  AFGL 
La  Junta.  Colorado  Selsmo-acoustlc  Monitoring  Array. 

Of  particular  Interest  for  the  vibro-acoustic  environment  of  the  site  is  the  fact  that  the 
array  lies  under  a  SAC  low-level  training  route  used  by  B-1,  B-52,  and  F-111  bomber  aircraft. 
In  this  route  SAC  aircraft  are  permitted  to  fly  as  low  as  130  meters.  The  La  Junta  site  Is 
essentially  In  the  middle  of  the  16.5  km-wlde  training  route  corridor.  However,  typical 
aircraft  flight  paths  take  them  within  approximately  1  km  east  or  west  of  the  site. 

3.3  Geologic  Setting 

The  site  Is  located  1300  m  above  sea  level  on  a  flat  plain  running  out  to  the  east  from  a 
SSW  trending  ridge  line  about  100  m  above  the  surrounding  terrain.  The  surface  layers  appear 
to  be  fine  grained,  sandy  soils  with  a  shallow  caliche  layer.  It  can  be  assumed  that  these 
layers  8ire  largely  the  deposits  from  the  erosion  of  the  ridge  line  to  the  west  of  the  site. 

Physlographlcally,  the  site  lies  on  the  western  margin  of  the  stable  continental  Interior  of 
North  America.*  Although  the  deeper  structure  has  not  been  studied  for  this  report.  It  can  be 
assumed  that  It  Is  typical  of  the  region,  consisting  of  largely  undisturbed  sedimentary  rocks 
laid  down  under  the  Cretaceous  epicontinental  sea.  Since  the  Devonian,  the  dominant 


*  Clark,  T.,  and  Steam,  C.  (1960)  The  Geological  Evolution  of  North  America,  The  Ronald 
Press  Co.,  New  York. 
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tectonic  process  for  the  region  has  been  a  general  uplift  associated  with  the  episodic  orgenles 
of  the  Rocky  Mountains. 


Figure  2.  Map  Showing  the  AFGL  Selsmo-acousUc  Array,  Local 
Topography,  and  Major  Cultural  Noise  Sources  In  the  Area. 


4,  GEOPHYSICAL  DATA  ACgUISITION  SYSTEM 


4.1  General  Description 


Tlie  selsmo-acoustlc  measurement  system  used  at  the  La  Junta  site  was  an  upgraded 
Standalone  Data  Acquisition  System  (SDAS)  as  described  by  von  Glahn.2  This  package,  the 


2  von  Glahn,  P.G.  (1980)  The  Air  Force  Geophysics  laboratory  Standalone  Data  Acquisiflon 
System:  A  Functional  Description,  AFGL  Report  No.  AFGL-TR-80-0317.  ADA160253. 


Geophysical  Data  Acquisition  System  (GDAS),  Is  a  unified  multi-channel  data  sampling, 
signal  conditioning,  and  r  ocesslng  system.  Data  used  in  this  report  were  acquired  with  the 
GDAS  recording  four  seismometers,  two  vertical  and  two  orthogonal  horizontal  Instruments, 
and  three  pressure  transducers.  All  of  the  sensors  were  collocated  for  estimation  of  system 
generated  noise.  (Collocated,  as  used  here.  Is  defined  to  mean  that  the  distances  betw'een 
sensors  is  less  than  0.1  times  the  shortest  wavelength  anticipated  within  the  system  pass 
band.) 

The  seismic  environment  was  measured  using  Hall-Sears  HS-IO-IB  velocity  transducers. 
These  Instruments  are  characterized  by  nominal  natural  frequencies  of  1  Hz  and  intrinsic 
sensitivities  of  the  order  of  1200  volts/ (meter/ sec).  Pressure  sensors  used  were  DJ  Instruments 
MLR  +/-  3450  Pascal  differential  transducers.  A  low  pass  anti-aliasing  Alter,  having  a  comer 
frequency  of  35  Hz  and  a  roll-off  of  36  dB/octive.  was  appAed  to  the  data.  This  provided 
nominally  flat  system  responses  between  1  and  30  Hz  for  both  seismic  and  pressure  signals. 

During  data  collection,  the  GDAS  was  programmed  to  record  30-sec  noise  samples  at  5- 
hour  intervals  for  all  available  channels.  A  sampling  rate  of  100  samples/sec  was  used  for  all 
channels  and  sample  sequences. 


4.2  Channel  Response  Calibrations 

The  seismic  channels  were  calibrated  in-situ  by  application  of  a  known  current  to  the 
calibration  colls  of  the  seismometers.  This  is  equivalent  to  the  application  of  a  known  force 
to  the  mass  of  the  Instrument.  Estimates  of  the  system  response  of  each  channel,  the  response 
due  to  the  instrument  and  signal  conditioning,  were  obtained  by  minimizing  the  least  squared 
error  between  the  observed  calibration  pulses  and  pulses  derived  from  theoretical  models  of 
the  system.  A  typical  response  curve  for  a  seismometer  is  shown  in  Figure  3  and  the  system 
parameters  for  each  seismic  charmel  are  given  In  Table  1 . 


Table  1.  Seismometer  Response  Characteristics 


Natural 

Intrinsic 

Channel 

Sense 

Frequency 

Damping 

Gain 

(Hz) 

(V/mm/sec) 

1 

Vertical 

1.096 

0.775 

172.5 

2 

Vertical 

1.100 

0.775 

201.7 

3 

N-S 

1.761 

0.381 

73.4 

4 


E-W 


0.650 


1.000 


509.9 
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Figure  3.  A  Typical  Seismometer  Channel  Response  Curve 
From  the  La  Junta  Selsmo-acoustlc  Array  (Channel  1). 


The  design  of  the  pressure  transducers  prevents  in-situ  calibration  of  the  complete  unit. 
Sensor  scale  factors  were  determined  by  loading  each  cell  with  a  pressure  defined  by  a  column 
of  water  of  known  height.  The  electronic  elements  of  the  channel  response  are  calibrated  In 
the  same  manner  as  the  seismic  channel  responses.  Figure  4  shows  a  typical  response  curve 
for  a  pressure  sensor  and  the  channel  response  parameters  are  given  in  Table  2. 


Table  2.  Pressure  Transducer  Response  Characteristics 


Natural 

Time 

Intrinsic 

Channel 

Frequency 

Constant 

Gain 

(Hz) 

(Sec) 

(V/mm/sec) 

5 

0.524 

1.7 

2.895 

6 

0.452 

1.7 

3.047 

7 


0.492 


1.7 


2.884 
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Figure  4.  A  Typical  Pressure  Sensor  Channel  Response  Curve 
from  the  La  Junta  Selsmo-acoustlc  Array  (Channel  7). 


4.3  System  Measurement  Errors 


Any  physically  realizable  measurement  system  will  record,  in  addition  to  the  desired 
signal,  some  level  of  noise.  This  noise  can  be  external  to  the  system  and  be  Introduced  as  an 
undesirable  signal.  In  addition,  the  measurement  system  will  Introduce  self-generated  noise 
Into  the  output  signal.  To  insure  the  fidelity  of  the  Input  and  recorded  signal,  the  system 
measurement  errors  must  be  significantly  below  the  amplitude  of  the  Input  signal. 

4.3.1  ADDITIVE  SYSTEM  NOISE 

An  estimate  of  the  additive  system  noise  Is  made  by  comparing  the  outputs  of  two  or  more 
collocated  sensors.  The  output  of  thejth  channel.  Sj(t),  Is  defined  as; 
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Sj(t)  =  Uj{t)  +  Nj(t) 


(1) 


where  Uj(t)  Is  that  part  of  the  signal  due  to  sensor  Input  and  Nj(t)  Is  that  part  due  to  system 
measurement  noise.  For  two  collocated  sensors.  J  and  k,  It  is  assumed  that  the  sensor  input 
term  is  common  to  both  output  channels  while  the  system  noise  term  is  Independent  and 
Incoherent  across  the  channels,  or: 


Uj(t)  =  u^(t) 


(2) 


and 


<Nj(t)  ■  Nk(t)  >  =  5jk 


(3) 


where  <  x  •  y  >  denotes  the  inner  product  of  x  and  y  and  is  the  Kronecker  delta  function. 
Then,  the  system  noise  spectrum  is  defined  as  the  output  signal  spectrum  minus  the  coherent 
spectrum  between  collocated  channels,  i  and  j.  Under  this  definition  a  lower  bound  estimate 
of  the  system  noise  will  be  made  as  any  coherent  noise  terms  are  assumed  to  be  part  of  the 
sensor  input  term. 

The  coherency  function  for  channels  1  and  J.  C,j(f).  is  given  by: 


C,j(f).|0,/f)j/jG„(f)-G„(f)p 


14) 


where  G,j  is  the  cross-spectral  density  function  between  channels  i  and  J  and  is  the  auto- 
spectral  density  function  for  channel  1.  In  turn,  the  Incoherent  noise  spectrum  is  given  by: 


G;(n=[i-c,^(f))  Gj,(t) 


IS) 


and  is  the  basis  for  the  estimate  of  the  GDAS  measurement  system  noise. 

Figure  5  shows  the  average  squared  coherency  function  for  the  collocated  vertical 
seismometers  (channels  1  and  2  at  the  La  Junta  site).  This  estimate  is  based  on  data  samples 
taken  from  all  67  available  time  windows.  Figure  6  shows  the  equivalent  function  for  two 
collocated  pressure  transducers  (channels  5  and  6).  A  generally  high  level  of  coherency  is  seen 
between  the  vertical  selsmlcs  in  the  band  of  approximately  15  to  35  Hz.  For  the  pressure 
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transducers,,  the  coherency  Is  much  lower  than  for  the  seismic  charmels  with  the  exception  of 
a  narrow  band  Just  above  10  Hz.  Figures  7  and  8  give  the  <ncoherent  spectral  estimates 
between  the  respective  sensors  based  on  the  67  data  segments.  These  plots  represent  the  best 
lower  bound  estimates  of  additive  system  generated  noise  for  seismic  and  pressure  charmels. 
The  additive  noise  is  approximately  4  to  6  orders  of  magnitude  in  power  above  the  theoretical 
digitization  noise,  shown  in  Figure  9,  of  a  15-blt  quantizer  (A/D),  as  used  in  the  GDAS.^ 
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Figure  5.  Average  Squared  Coherency  Function 
for  Two  Collocated  Vertical  Seismometers  from  the 
La  Junta  Selsmo-acoustlc  Array  (Charmels  1  and  2). 


®  Bermett,  W.  (1948)  Spectra  of  Quantized  Signals,  Bell  Systems  Technical  Journal  27 
(No.  3):446-472. 
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Figure  6.  Average  Squared  Coherency  Function  for 
Two  Collocated  Pressure  Trauisducers  from  the 
La  Junta  Selsmo-acoustic  Array  (Channels  5  and  6). 
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Figure  7.  Average  Incoherent  Spectrum  for 
Two  Collocated  Vertical  Seismometers  at  the 
La  Junta  Selsmo-acoustlc  Array  (Channels  1  and  2) 
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Figure  9.  Theoretical  Digitization  Noise  of  a  15-blt 
Quantizer  Representing  the  GDAS  A/D  System. 


4.3.2  SIGNAL  INDUCED  NOISE 


One  potentially  significant  source  of  signal  induced  noise  is  the  folding  or  aliasing  of  out- 
of-band  power  back  into  the  measurement  data  band.  The  ability  of  the  measurement  system 
to  mitigate  this  effect  is  characterized  by  the  protection  ratio,  the  ratio  of  the  in-band  system 
response  to  the  response  at  the  first  fold-aliased  frequency.^  This  can  be  expressed  as: 


P,(f)-T,(f)/T,(t3-t)  0<f<l„ 


where  T,(0  is  the  system  response  of  channel  1.  fs  is  the  sampling  frequency  and  f^  is  the 
Nyqulst  fi-equency  or  O.Sfg.  The  protection  ratio  for  chamiel  1,  a  vertical  seismometer,  is 
shown  in  Figure  10  and  for  channel  7,  a  pressure  transducer,  in  Figure  11.  Given  equal 
amplitude  inputs  to  the  measurement  system  at  some  in-band  frequency  and  its  first  fold 


*  Blackman.  R.,  and  Tukey,  W.  (1958)  The  Measurement  of  Power  Spectra,  Dover  Publication 
New  York 


Figure  1 1 .  Protection  Ratio  for  Channel  7,  a  Pressure 
Transducer  of  the  La  Junta  Selsmo-acousUc  Array. 


4.4  Data  Quality 


The  quality  of  the  data  recorded  at  the  La  Junta  array  can  be  quantified  In  terms  of  the 
signed  to  noise  ratio  of  the  recordings.  There  are  two  potentially  significant  noise  sources 
affecting  the  data,  the  additive  systems  noise  and  out-of-band  signals.  It  is  desirable  to 
consider  two  forms  of  the  signal  to  noise  ratio  to  isolate  the  contribution  of  each  noise  source. 

In  Section  4.3.1,  the  signal  was  defined  to  be  coherent  across  collocated  sensors  and  the 
additive  system  noise  was  defined  as  the  Incoherent  signal.  Thus,  the  signal  to  noise  ratio 
assuming  only  additive  system  noise  is  given  by  the  ratio  of  the  coherent  spectra  to  the 
Incoherent  spectra  for  any  given  channel,  or: 


S/N,=[G;(f)/(^(f)j2 


(7) 


or 
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PROTECTION  RATIO 


S/N,  =  [Cy(f)/{l-Cy(f)}j2 


(8) 


for  channel  i.  Figures  12  and  13  show  the  signal  to  noise  ratios  for  vertical  seisniics  and 
pressure  transducers,  respectively. 
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Figure  12.  Signal  to  Noise  Ratio  (S/N) 
Assuming  Only  Additive  System  Noise 
for  a  Typical  Seismic  Channel  Recording. 
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SIGNAL  TO  NOISE  RATIO 


FREQUENCY  (HZ) 

Figure  13.  Signal  to  Noise  Ratio  (S/N)  Assuming 
Only  Additive  System  Noise  for  a  Typical  Pressure 
Channel  Recording. 

Alternatively,  the  signal  to  noise  ratio  can  be  described  as  the  ratio  of  the  maudmum  In- 
band  signal  amplitude  to  the  additive  measurement  noise  plus  the  first  fold  aliasing  error. 
This  version  of  the  signal  to  noise  ratio  Is  given  by: 


S/N',  .  [c;(f)/{G",(f)  H.  G’(f„)  pf  (0}]^ 


where  G*',(fN)  Is  the  maximum  coherent  signal  power  at  the  Nyqulst  frequency,  ffj.  No 
significant  differences  were  noted  between  the  two  definitions  of  signal  to  noise.  This 
Indicates  that  the  Inband  signals  are  not  being  distorted  by  aliasing  of  power  from  outside  the 
measurement  band. 
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5.  SEISMO-ACOUSTIC  NOISE  ENVIRONMENT 


5. 1  Introduction 

Seismic  noise  at  any  site,  particularly  for  shallow  Instruments,  can  generally  be 
characterized  as  a  combination  of  weak  propagating  disturbances  from  distant  and  near 
sources  often  associated  with  cultural  activity  and  meteorological  conditions.  Pressure  noise 
results  from  the  pressure  fluctuations  due  to  turbulent  motions  within  the  atmosphere,  that  Is, 
wind  generated  noise.  It  Is  well  known  that  these  noise  fields  demonstrate  strong  seasonal 
and  dally  variation  due  piim.-irlly  to  meteorological  effects.®-  ®  It  was  the  goal  of  this  study  to 
evaluate  parameters  to  describe  the  statistical  character  of  the  noise  fields  observed  at  the  La 
Junta  array  and  to  establish  the  degree  of  correlation  between  the  pressure  and  seismic  noise 
signals. 

5.2  Broadband  Characteristics 

Dur^iig  the  period  9  through  29  October  1986,  the  La  Junta  selsmo-acousUc  array  was 
triggered  every  5  hours  to  obtain  30-sec  segments  of  data  at  100  samples/sec.  A  total  of  74 
data  segments  were  collected  during  this  window.  The  raw  data  files  were  reviewed  and  5.12- 
sec  samples  selected  for  processing  to  eliminate  data  either  demonstrating  equipment 
malfunctions  or  containing  apparent  discrete  transient  events.  Of  the  original  73  data 
segments,  67  segments  were  used  for  analysis.  Figure  14  shows  typical  examples  of  the  seismic 
and  pressure  noise  data  as  recorded  at  the  La  Junta  array. 

Overall  standard  deviations  for  each  channel  were  2.94  x  10”®  and  8.61  x  10”®  mm/ sec  for 
the  vertical  and  horizontal  seismic  components,  respectively,  and  1.96  x  10“^  Pa  for  the 
atmospheric  pressure  noise.  Figure  15  displays  the  cumulative  probability  distribution 
functions,  plotted  on  normal  probability  scales,  of  the  vertical  and  horizontal  seismic  and 
pressure  noise  amplitudes.  The  strong  divergence  from  linear  behavior  In  the  tails  of  the 
distributions  demonstrates  the  non-Gausslan  nature  of  the  broadband  seismic  and  pressure 
noise  fields.  The  non-Gausslan  behavior  of  the  data  Is  attributed  to  non-statlonarity  of  the 
distribution  parameters. 


®  Akl,  K.,  and  Richards,  P.  (1980)  Quantitative  Seismology:  Theory  and  Methods.  W.  H.  Freeman 
and  Company,  San  Francisco. 

®  Rlngdal,  F.,  and  Bungum,  H.  (1977)  Noise  level  variation  at  NORSAR  and  Its  effect  on 
detectability.  Bull  Seism.  Soc.  Amer.,  67(No.  2):479-492. 
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Figure  15a. 
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Figure  15.  Cummulatlve  Probability  Distributions  for  Vertical 
(a),  and  Horizonteil  (b)  Seismic  and  Pressure  (c)  Noise  Amplitude 
Data  on  a  Normal  Probability  Scale. 


In  addition  to  the  Gaussian  distribution,  the  data  was  compared  to  the  Gumbel  Type  II 
extreme-value  distribution.  Distribution  parameters  were  obtained  based  on  the  observed 
maximum  absolute  value  of  each  measured  field  in  1.01 -sec  windows.  The  parameters  of  these 
distributions  are  (u  =  3.8  x  10”®,  K  =  2. 56]  and  |u  =  9.1  x  10”®,  K  =  2.29)  for  vertical  and 
horizontal  selsmics,  respectively,  and  Ju  =  0.108,  K  =  2.12]  for  the  atmospheric  pressure  field. 
The  observed  and  theoretically  estimated  distribution  functions  are  shown  In  Figure  16. 

While  the  estimated  functions  are  not  perfect  representations  of  the  observed  distributions, 
vertical  seismic  and  pressure  distributions  are  considered  to  be  adequately  fit. 


22 


I  T  t  T  I  I  I  I  I  [  I  I  I  I  I  I  I  I  I  I  T  I  -i-ri  ri  11  I  1  T  I 

0.0001  0.0002  0.0003 

AMPLITUDE  IN  MM/SEC 


TrjXTTTT 

0.0004 


Figure  16a. 


r 

0.0005 


23 


AMPLITUDE  IN  PA 


Figure  16c. 

Figure  16.  Observed  (SoUd)  and  Fit  (Dashed)  Gumbel  Type  II  Extreme 
Value  Distributions  for  (a)  the  Vertical  Seismic  Noise,  (b)  Horizontal 
Seismic  Noise,  and  (c)  the  Atmospheric  F>ressure  Signal. 


Figures  17  through  19  summarize  the  temporal  variation  of  maximum  magnitude,  the 
largest  absolute  value  of  the  signal  amplitude  and  the  sample  standard  deviation  for  the 
vertical  and  horizontal  components  of  seismic  motion  and  a  pressure  channel  during  the 
observation  period.  Each  vertical  har  on  the  plots  represents  the  value  obtained  from  one 
5.12-sec  data  sample  with  a  5-hour  time  shift  between  sequential  bars.  The  cross-hatched  bars 
fill  for  samples  that  were  eliminated  from  analysis  due  to  apparent  equipment  failures  or 
obvious  transient  events. 
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Figure  17.  Amplitude  Statistics  of  the  Vertical 
Seismic  Noise  Field  at  5-Hour  Intervals. 
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Figure  18.  Amplitude  Statistics  of  the  Horizontal 
Seismic  Noise  Field  at  5- Hour  Intervals. 


26 


MM /SEC  MM/SEC 


ABSOLUTE  MAXIMUM 


Ihui  -  _ _ 

” 

s 

's 

s- 

s 

s 

3 

s 

s 

N 

'N 

V 

s 

S 

S 

V 

s 

s 

N 

•v 

S. 

V 

Ns 

s 

S 

N 

s 

S 

S 

N 

N 

1  !  1 - 1 - 1 - 1 - 1 - 1 - 1 - i  i - 1 - 1 - 1 - 1 

- 1  i - 1 - 1 - 1 - ^ - i - 1 - 1 - 1 - 1 - i - 

1  1 — \ — 1 

4. 


0.03 


STANDARD  DEVIATION 


0.  110.  220.  330. 

TIME  IN  HOURS  -  (OFFSET  FROM  START  TIME) 


Figure  19.  Amplitude  Statistics  of  the 
Pressure  Noise  Field  at  5-Hour  Intervals. 


Considering  each  plot  separately,  It  Is  seen  that  the  maximum  magnitude  and  the  standard 
deviation  on  each  plot  move  essentially  In  lock  step.  In  other  words,  Increases  In  maximum 
magnitude  are  associated  with  Increases  In  the  general  level  of  the  noise  and  are  not  the  result 
of  the  injection  of  Isolated  spikes  within  an  essentially  constant  background. 

Contrasting  the  figures.  It  Is  apparent  that  there  is  some  degree  of  correlation  between  the 
pressure  noise  and  surflcial  noise  levels  and  that  this  correlation  is  most  prominent  during 
the  period  between  30  emd  60  hours  Into  the  sampling  window.  The  strength  of  the  correlation 
during  this  period  Is  emphasized  due  to  the  extremely  high  pressi>re  noise  level  observed 
which.  In  turn,  de-emphaslzes  the  variation  in  pressure  noise  levels  at  other  times.  While  a 
careful  examination  of  the  data  outside  this  window  shows  correlation  between  periods  of 
high  seismic  and  pressure  noise.  It  Is  also  apparent  that  It  Is  not  strong,  similarly,  a  general 
correlation  can  be  seen  between  the  vertical  and  horizontal  component  of  the  seismic  noise. 

As  Is  often  found  for  seismic  noise,  however,  the  vertical  component  demonstrates  a  smaller 
degree  of  variation  than  do  the  horizontal  motions.® 

Figure  20  Is  a  polar  plot  of  the  relative  frequency  of  the  azimuth  of  particle  motion  in  the 
horizontal  plane  based  on  azimuth  estimates  from  all  67  data  segments  In  the  analysis.  It  is 
apparent  from  this  plot  that  the  horizontal  particle  motion  distribution  is  non-isotroplc  and 
is  biased  along  an  axis  between  NE-SW  and  E-W,  approximately  80°-260°.  No  further  attempt 
was  made  to  Identify  the  cause  or  In  any  other  way  examine  this  feature  of  the  seismic  noise 
field. 
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Figure  20.  Relative  Frequency  of  Horizontal 
Particle  Motions  Recorded  at  the  La  Junta  Site. 

5.3  Spectral  Characterization 

Mean  periodograms  evaluated  for  vertical  and  horizontal  selsmlcs  and  pressure  noise  data 
using  the  ensemble  of  sample  segments  are  shown  In  Figure  21.  The  pressure  perlodogram 
shows  a  typical,  nearly  linear,  wind  noise  roll  off  out  to  the  anti-aliasing  filter  cut-off 
frequency  of  30  Hz.^  At  low  frequencies,  below  approximately  1  Hz,  the  seismic  periodograms 
show  significant  power  commonly  associated  with  the  storm  mlcroselsm. 


^  Bruce.  R.  (1971)  Field  measurements:  equipment  and  techniques.  In  Noise  and  Vibration 
Controi  L.  Beranek,  ed.,  McGraw-Hill  Book  Company.  New  York. 
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Figure  21.  Mean  Perlodograms  of  the  Vertical  (a)  and 
Horizontal  (b)  Seismic  and  Pressure  (c)  Noise  Data. 

As  was  demonstrated  In  Section  5.2,  the  broadband  seismic  and  pressure  noise  field  cannot 
be  described  as  Gaussian.  However,  that  does  not  preclude  the  possibility  that  within  specific 
spectral  bands  the  fields  are  Gaussian.  The  hypothesis  that  the  noise  fields,  within  narrow 
bands.  Is  Gaussian  can  be  tested  using  a  technique  proposed  by  Sax.®  This  test  is  based  on  the 
fact  that  periodogram  coefficients  for  a  stationary,  Gaussian  process  are  distributed  as  a 
Rayleigh  process.  Then,  If  the  observed  coefficients  are  plotted  against  a  Rayleigh 
distribution,  a  straight  line  fit  of  the  coefficients  should  be  observed. 

Figure  22  shows  examples  of  the  test  results  for  the  vertical  and  horizontal  selsmlcs  and 
the  pressure  data  for  a  band  centered  at  a  fi-equency  of  5.1  Hz.  In  this  case,  as  for  all  other 
frequencies  tested,  the  results  are  distinctly  non-linear  and  do  not  support  an  assumption  of 
an  active  Gaussian  process. 


®  Sax.  R.  (1968)  Statlonarity  of  seismic  noise.  Geophysics,  33  (No.  4):668-674. 
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Figure  22.  Distribution  of  Periodogram  Coefficients  at 
5.1  Hz  Versus  a  Rayleigh  Distribution  for  Vertical  (a) 
and  Horizontal  (b)  ^ismlc,  and  Pressure  (c)  Noise. 


5.4  CrosS'Fleld  Coherency 


Figures  23  through  25  show  the  coherency  squared  function  across  noise  fields,  specifically 
between  vertical  and  horizontal  seismlcs,  vertical  selsmlcs  and  pressure  and  horizontal 
selsmlcs  and  pressure.  In  contrast  to  the  like  field  coherencies,  such  as  shown  in  Figure  5  for 
two  vertical  seismometers,  only  very  weak  coherencies,  at  best,  are  found  across  fields. 
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Figure  25.  Coherency  Squared  Function  Between 
Horizontal  Seismic  and  Pressure  Noise. 


6.0  CONCLUSIONS 


Analysis  of  seismic  and  atmospheric  pressure  noise  environment  data  taken  during  the 
period  9  to  29  October  1986  at  the  La  Junta.  Colorado  array  site  demonstrates  that  p>eak 
ambient  conditions  for  this  site  and  during  this  window,  needed  to  define  event  detection 
capabilities,  can  adequately  be  modeled  as  Typ>e  II  Gumbel  extreme-value  distributions. 
Further,  It  has  been  shown  that  there  is  at  best  only  weak  correlation  between  the  various 
noise  fields  studied  in  both  a  broad  and  narrow-band  sense. 
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